Iranian Journal of Numerical Analysis and Optimization 


Vol. 13, No. 4, 2023, pp 747-762 — 
https: //doi.org/10.22067 /ijnao.2023.78096.1165 _ 
https://ijnao.um.ac.ir/ 


How to cite this article ro} 
BY NC Research Article 


Optimal control analysis for modeling 
HIV transmission 


K. R. Cheneke 


Abstract 


In this study, a modified model of HIV with therapeutic and preventive 
controls is developed. Moreover, a simple evaluation of the optimal control 
problem is investigated. We construct the Hamiltonian function by way 
of integrating Pontryagin’s maximal principle to achieve the point-wise 
optimal solution. The effects obtained from the version analysis strengthen 
public health education to a conscious population, PrEP for early activation 
of HIV infection prevention, and early treatment with artwork for safe 
life after HIV infection. Moreover, numerical simulations are done using 
the MATLAB platform to illustrate the qualitative conduct of the HIV 
infection. In the end, we receive that adhering to ART protective prone 
people, the usage of PrEP along with different prevention control is safer 
control measures. 


AMS subject classifications (2020): Primary 45D05; Secondary 42C10, 65G99. 
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1 Introduction 


Human immunodeficiency virus (HIV), the cause of HIV infection, has no 
curative medication until now [2]. Moreover, the long-time existence of the 
virus in the body leads to a serious infection called acquired immunodefi- 
ciency syndrome (AIDS) disease [6]. However, optimal controls are the ef- 
fective way to combat HIV transmission and progression in the community 
(2, 6]. Public health education, condom, and anti-retrovirus therapy are the 


Received 5 August 2022; revised 26 May 2023; accepted 16 June 2023 


Kumama Regassa Cheneke 
Department of Mathematics, Wollega University, Nekemte, Ethiopia. e-mail: kuma- 
maregassa@gmail.com 


TAT 


Cheneke 748 


major measures taken by both governmental and nongovernmental institu- 
tions to stop further progression and transmission of HIV in the populations 
[1, 3, 12, 4, 5, 7, 8, 9]. Moreover, effective pre-exposure prophylaxis (PrEP) 
is the drug used to prevent the survival of HIV in the human blood [14]. 
On the other hand, the abstinence of sexual practices through the activation 
of public health education reduces the fate of acquiring HIV from poten- 
tially infectious individuals. Mathematical models are very important tools 
to describe the behavior of biological events. Particularly, with the great 
contribution of Pontryagin’s Maximum Principle (PMP) in the construction 
of optimal control problems, the nature of biological dynamics is studied in- 
tensively [13, 14, 15, 18, 19, 20, 10, 22, 23]. Based on the works done in [12], 
the motivation of this study is due to the significant contribution of pub- 
lic health education and prophylaxis in controlling the transmission of HIV 
infection among human individuals. Particularly, abstinence due to consol- 
idated public health education builds positive awareness toward controlling 
oneself, whereas prophylaxis helps to prevent the progression of HIV in the 
human body. Mathematical models are important tools to control infections 
(24, 11, 27, 26, 16, 17]. In this study, we have included prophylaxis, antiretro- 
viral therapy (ART), and prevention for controlling the transmission of HIV 
infection by modifying the model studied in [12]. 


2 Formulation of model 


In this study, a mathematical model is formulated by classifying the total 
population into compartments of (i) Susceptible individuals (S), (ii) Individ- 
uals on Pre-exposure prophylaxis (E), (iii) HIV infected with primary stage 
(P), (iv) Not on treatment HIV infected individuals (J), (v) HIV Undetectable 
individuals (U), and (vi) On treatment HIV infected individuals (1). 

Moreover, the subsequent assumptions are considered in the modeling of 
the infection (i) a new susceptible individuals becomes susceptible at recruit- 
ment rate of A, (ii) individuals in S transfer to E due to taking PrEP at the 
rate of p; (iii) transmission rate of HIV infection from individuals in P to S is 
(1 and transmission rate of HIV infection from I to S is 69; (iv) individuals 
transfer from P to I at progression rate of €; (v) individuals transfer from P 
to J at transfer rate of 7; (vi) individuals transfer from J to I at transfer rate 
of y; (vii) individuals in the compartment J die due to infection at the rate 
¢; (viii) individuals transfer from I to U due to adherence to ART at trans- 
ferring rate of 0; (ix) individuals transfer from U to I due to default using of 
ART at the rate of ¢; (x) natural induced death rate of all people is ju; (xi) 
AIDS induced death rate is 6; (xii) in this study, standard incidence rate is 
applied; (xiii) PrEP engagement effort is u3; (xiv) Condom using effort is u1; 
(xv) ART using effort is ue. 

The pictorial representation of the deterministic model with control mea- 
sures is given in Figure 1. 
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Figure 1: Schematic diagram of HIV transmission dynamics. 


The deterministic model of population dynamics subject to HIV infection 
in the presence of control measures is given by 


= See (usp 28, 
dE 
He 7 UPS — (ut eye 
2 = i ak 
(1) 
dJ 
7 ie Ce a oLe 
SHEP HS + OU ~ (uO tu +6)I 
= wz6I ~ (9+ u)U 


with initial conditions: $(0) > 0, 
t= 1,2, 


E(0) > 0, P(0) > 0, J(0) > 0, 1(0) > 
0, U(0) >0, 0< uy <1, 4. 


3, 


3 Analysis of the model without control 
3.1 Invariant region 


Theorem 1. The solution of model (1) is invariant in the region 2 proper- 
subset of six-dimensional space over the set of nonnegative real numbers such 
that 


= {(S,E, P, J,U,1) € R§.: N(0) < 


Ely 


}. (2) 


Proof. The equations of model (1) gives the subsequent equation: 
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dN 
=X N— ol J. 
a a CJ, 
which implies 
dN 
— <A—-wpN. 
da-° » 


Applying mathematical procedures, the preceding inequality gives 
r aN 
N(t)< —- (2 — n(o) em, 
wo \E 


which implies, as time ¢ varies, the total population size is bounded for all 
time t, with the given initial condition. 


3.2 Nonnegative property 


Theorem 2. All solution variables of model (1) without control are nonneg- 
ative in the stated invariant region of the solution. 


Proof. Consider the first equation of model (1) without control. Then 


dS S(Bi1P + Bol + B37) 
dt N 


t PE — pS, (3) 
which implies 
dS _S(6iP + Bal + Bs) 
dt ~ N 
Solving the preceding inequality, we get 


Ls. (4) 


t (81 P(€)+621(€)+83 J(€)) 
pot— fy RE} 2 


S(t) > S(O)e (5) 


Hence, based on the initial condition, the susceptible population size is non- 
negative for all time t. 


3.3 Basic reproduction number 


The basic reproduction number Ro of model (1) without control is the average 
number of infected individuals produced by typical infectious individuals in 
the susceptible population during the entire period of infection. Based on the 
techniques applied, we compute basic reproduction numbers from model (1) 
without control as follows. Let F and V be the Jacobian matrices obtained 
from model (1) as given below: 
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The spectral radius p (F V-t) computed from next-generation matrix 
FV~! of foregoing matrices is given by 


= By i Bon 
chy (Pew ele) 
Paleur ey EC Hy) 
(6+uE+n+m(yt+o+u) 


p(FV~") 


Therefore, by the definition, we obtain 


- fon Bon 
Cty Mage wre Cay) 
B3x(Eu+ Ey + &C +07) 

(O+my(E+tn+u)(y+o+p) 


Ro 


3.4 Global stability of disease-free equilibrium 


Theorem 3. The global stability of a disease-free equilibrium point is de- 
scribed as a steady state where the trajectory of solution shows the tendency 
of moving toward it for all time tf. 


Proof. To show the global stability of disease-free equilibrium, we incorporate 
the method applied in the works of [21]. Next, from the computed matrices 
for construction of next-generation, we obtain 


By Bo B3 0 aaa 0 0 0 
{0000 _ —yn ytCtp 0 0 
a 0000)’ oo —€ —y éd+p —-¢ 
0000 0 0 0 d+p 


Moreover, the rate of change of variables (P, J,I,U) at disease-free equilib- 
rium can be written as 


aE P 
as J 
dts 

Tr U 


Therefore, by the comparison method applied in [21], we justify that 
model (1) without control has a globally asymptotically stable disease-free 
equilibrium. 
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4 Extension to control problem 


The deterministic model of population dynamics subject to HIV infection, in 
the presence of control measures, is given by 


8 a ieee ey) pm bape ihe. 
dE 
He UsPS — (H+ v) EB, 
“= Sse (E+n+n)P, 
(6) 
dJ 
ge TOP +S + wd, 
OS EPL THOU (29 +54 WL, 
dU 
ai = u6l — (d+ p)U, 


with initial conditions $(0) > 0, E(0) > 0, P(0) > 0, J(0) > 0, I(0) > 
0, U(0) >0, 0< uj <1, 1=1,2,3. 

To study the optimal levels of the controls, we define the Lebesgue measurable 
control set U as 


U = {(uy, v2, 3) :0< uy < 1,0 < ug < 1,0 < ug < LO< tate (7) 


Our goal is to find the optimal controls that minimize objective functional J 
given by 


tf 

J= min | oP +col +co3J + (wiuz + wous + w3uz), — (8) 
(u1,u2,u3) Jo 2 

where cj, j = 1,2,3 and w; ,2 = 1,2,3 are constants. The expressions 

0.5wju?,i = 1,2,3 are costs associated with controls. The form of cost is 

quadratic because we assumed it to be nonlinear in nature [24]. Also, for 

four optimal controls uj, u3,u3, we have 


JI(ujz,u5,u3) = min{ J(u, ua, ug) : U1, U2, u3 € US, 


where U = {(u1, U2, us): 0 < ur < 1,0 < ue < 1,0 < ug < 1}. Furthermore, 
uy ,, uzandu3 are measurable controls. 
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4.1 Existence of optimal control solution 


Theorem 4. The optimal control solution of a control problem exists if the 
following Fleming’s and Rishel’s conditions are satisfied: 


(i) The set of all solutions to optimal control problem and objective func- 
tional must be nonempty. 


(ii) The state system is a linear function of controls with coefficients de- 
pendent on state variables and time. 


(iii) The integrand in objective functional is convex and bounded above by 


dy (|uy|?+]u2|?-+|u3|7)¢—de < cyP+c2I+e3I+4 (wyut + woz + wu) ,d, > 


Oandd>1. 


Proof. We employ the method from to demonstrate the existence of opti- 
mal control. The condition (i) is satisfied if the state system has bounded 
coefficients. Additionally, the state system operates in accordance with 
controls, satisfying requirement (ii). The integrand in the objective func- 
tional is used to demonstrate condition (iii). Moreover, c1P + col + c3J + 
1 


7 (wi us + woud + w3u3) is convex on U as any constant, linear and quadratic 


are convex. Furthermore, assume that there are d,,d2 > 0, andd > 1 satisfy- 


ing dy(|ui|?+|u2|?+|u3|7)?—de < cy P+cgI+c3I+4 (wyut + wouz + w3us) ,dy = 


min {w;,i = 1,2,3}, d = 2, and dz is the half of coefficient of control func- 
tions. Therefore, the optimal solution exists. 


4.2 The Hamiltonian and optimality system 


The PMP stated the necessary conditions that are satisfied optimal pair. 
Hence, by this principle, we obtain the Hamiltonian function (H) defined as 
[24] 


1 
A(a2,u,t) =a P+ cI +e3J + 3 (wiuz -f wus + w3u3) 


dI |) au 
dt | °° dt’ 


dP 
+ A3 t +4 


dJ 
7 r +As5 


dS dE 
+ r1 a + r2 Fi 


dt 
where ;,2 = 1,2,3,4,5,6 are the adjoint variable corresponding to state 
variables S, FE, P,J,I, and U, respectively, and to be determined using the 
PMP for the existence of optimal pairs. 


Theorem 5. Let S,F,P,J,I,U be optimal state variables and let optimal 
control u;,2 = 1,2,3 be the optimal controls. Then there exist costate vari- 
ables M1, r2, A3, a; As, andAg that satisfy 
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dA, OH dr2 OH 


dt Os’ dt OE’ 
d\3 OH d\y ‘OH 
dt OP’ dt =OJ’ 
d\s _- OH d\y OH 
dt Ol’ dt OU’ 


with transversality or final time conditions A; (tr) = A(t) = As(tp) = 
Aa(tp) = As(t¢) = Ao(t¢) = 0, where H is the Hamiltonian function. More- 
over, the optimal controls uj, u4, and ug are uj = min{0, max{ 25/02—M) 1}} 


iN) ? 
and us = min{0,max{ As) 1}}, over the constraints 0 < u, < 1,0 < 


ug <1. 


Proof. The PMP gives the standard form of adjoint equation with transver- 
sality conditions. Now, differentiating the Hamiltonian function with respect 
to state variables, we have 


dd oH BP + Bol + B3J)N — S(BiP + Bol + Bad 

7 = ae =(1 ie * a) Ne 2 2 3 Oy A3) 
+ U3/P (Aq = A2) + LA, 

dX OH 

a ap = (A2 — A1) + Aa, 

d\s 8H (1 —t) (BSN — 8 (fP + Bol + Bal) 

at = 5p Se U1 1 a 2 3 (255) 
+ €(A3 — As) +9 (Az — Aa) + HAs, 

d\y_ OH (1—wu1)(63SN —S(fiP + Bol + B3J)) cme 

dt Os N2 oe ee 

(Aa —A5) + (C+ BAS, 

dX OH 1-— BoSN — S(6,P + Bol + b3J 

“ --3 re ur) (Be a 2 3 Yor, ds) 
+ u20(A5 — Ag) + (6+ LAs, 


dy OH 


“a ou (Ae — As) + HAs. 


Furthermore, the characterization of optimal controls uj,us and u3 shows 
that 


OH OH OH 


Ou Oug Ou3 =o 


Hence, optimal controls over 0 < uy < 1,0 < ug < 1,0 < ug < 1 are given by 


: 5S (81P + Bol + B3J) (Ag — A1) 


Uy =U = 
1 ‘l 
w,N : 
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en SI OGEIS) 
ee pg 
»_, _ PS (Ai — Az) 


Therefore, the bounds of the optimal control variables are given by 


SP thatthe DOs) gg < SrPtPaltPeDOs—M) < 1, 
ur =< 0 wiN if S(O: P+ 821-8 )09= M) <9, 
: ifl< SOP Hat Bad) Ax) 


AI(As—r6) if Q Oe) 2, 
us=20 if uO Aw) <0, 
1 ifl< G1(As—do) 

W2 


} 


S(Ai—A2) S(A1—Az) 
oo 


us = 4 0 if Bra) < Q, 
1 if 1 < 250=N) a 


In a compact form, the optimal controls can be written as 
ut = min{0, max{ S(LP+Bal+ Bs DOs—M) 14 


us = min{0, max{ 1As—As) 1h}, 


U3 = min{0, max{ 251-2) 1}}. 


Moreover, the optimality system of the optimal control problem can be 
written as 


dS l-—wu S(6,P + Pol + b3J 
t =i ( v ( hi Z . WE (usp t L) S, 
d 


B 
qe u3PS — (ute) B 
dP (1-m)S(AiP + Bol + B3J) 


7 N (Erg ee, 

dJ 

oe (aah hp), 

EP +4 +0 ~ (ws +544) T 

dry «gy EE SN ee Pal a) 


=(1 
dt ( N? 


+ U3/P (Ar = A2) + LA, 


(A1 — A3) 
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“ = (Ao — A1) + Ad, 

“ _ C1 rf (1 = U1) (3,SN = SAP + Bol + B3J)) Oy de) 
+ €(A3 — As) +7 (Az — Aa) + HAs, 

dry _ (1-11) (B3SN — S(GiP + Bal + B3J)) 

dt —_ N2 (At A3) 

(Aa — As) + (C+ B)Ad, 

os = Gs + (1 = u1)(B2SN = SPs Bol + B3J)) 1 es 
ae u29(A5 = Ae) + (0 + 1)A5, 

dX 

a =$(As — As) + HAs, 


with Ai (ts) = Ag(te) = As(te) = Aa(ty) = Aste) = Ao(te) = 0,S(0) = 
So, P(0) = Po, J(0) = Jo, 1(0) = Ip,U(0) = Uo. 


4.3 Numerical simulations and discussion 


4.3.1 Analysis using the numerical methods 


In this study, the numerical methods are involved in simulating the general 
results of the analytical findings that give real meaning to both the mathe- 
matical and biological communities. Furthermore, the parameter values used 
in the simulation are either taken from the literature or assumed, as given 
in Table 1. Also, wy, = 50, we = 20, w3 = 30,c1 = 5, co = 25,T = 20, 5(0) = 
1000, H(0) = 0,W (0) = 300, /(0) = 500, U(0) = 0, A(0) = 0. 

Moreover, MATLAB software is applied in the simulation process. Fractional 
derivatives and stochastic findings are widely applied as reviewed in this pa- 
per. Hence, we incorporate both forward and backward sweep methods of 
fourth-order Runge-Kutta method to simulate the results. The applied con- 
trol strategies are as follows: 

Strategy 1: Using together control measures u; and ug. 

Strategy 2: Using together control measures u; and u3. 

Strategy 3: Using together control measures ug and u3. 

Strategy 4: Using together control measures uj, u2, and us. 

Moreover, we have used the parameters given in Table 1 to simulate subse- 
quent numerical solutions. 
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Table 1: Parameter/constants value. 


Parameter/constants Value 


200 
Bi 0.9915 
Bo 0.75 
Bs 0.9815 
E 0.5 
LU 0.02 
n 0.5 
¢ 0.1 
o) 0.09 
0 0.5 
é6 1 

p 0.1 
Y 0.1 
~ 0.001 


Based on the aforementioned control strategies, the following numerical 
simulations are performed. 


t(years) 


Figure 2: Primary HIV infected population. 
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Figure 3: HIV not tested population. 


=0,u,=0 


2 #0,u, #0 


0 5 10 15 20 
t(years) 


Figure 4: HIV infected and on treatment individuals. 
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Figure 5: Control functions effect illustration. 
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Figure 6: Adjoint variables condition descriptions. 


4.3.2 Numerical results and discussion 


This study develops and analyzes a mathematical model of HIV with the best 
possible control measures. The conceptual diagram for population dynamics 
is shown in Figure 1. The inclusion of the intervention with control serves 
to emphasize the significance of control measures in minimizing the effects of 
HIV infection. The numerical simulation results are shown in Figure 2 and 
show how a successful combination of control methods lowers the number 
of newly infected people. The numerical results in Figure 3 show that a 
reduction in the number of people who have not begun ART is shown when 
an intervention with control functions is present. Figure 4 shows a simulation 
of the number of HIV-positive people who are now receiving treatment. The 
results show that the intervention with three control groups dramatically 
lowers the number of people infected with HIV. When applied correctly, uj; 
and ug are effective from the beginning to the end of initiation, as seen in 
Figure 5, where the applied control functions are simulated. Control ug, on 
the other hand, makes a smaller contribution to regulating HIV infection 
dynamics because of its limited availability. In Figure 6, the adjoint variable 
is simulated to show that the transversal requirement has been satisfied. 


5 Conclusion 


According to the results of analytical and numerical simulations, adopting 
the best control measures to stop the further progression and transmission 
of HIV dynamics is more successful if done before the HIV infection even 
begins to spread. Additionally, maintaining ART and protecting those who 
are susceptible are considered the most crucial ways to lessen the effects of 
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HIV infection. Intervention with pre-exposure prophylaxis contributes less 
to lowering the risk of HIV infection since it is less affordable and accessible. 
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